MODELING AN ELECTRONIC DEVICE 

BACKGROUND OF THE INVENTION 

Modeling of electronic devices e.g. for simulation of electronic circuits 
comprising such devices, as known for example from the SPICE circuit 
5 simulation system, requires to know or reveal details about the internal 
construction of the device to be modeled. Often manufacturers of devices (like 
operational amplifiers) provide model parameter to be used in the S PICE 
system. 

SUMMARY OF THE INVENTION 

10 It is an object of the invention to provide an improved modeling of electronic 
devices. . The object is solved by the independent claims. 
Preferred embodiments are given by the dependent claims. 

According to the present invention it is not necessary to provide an analytical 
model describing the internal construction, but the measured (external) 
15 behavior of the device in time or frequency domain is used. In a preferred 
embodiment, the measured electrical behavior comprises a measured signal 
response on a predetermined electrical signal input to the device, and 
preferably comprises a measured signal response on a step signal. 

The device can contain any kind of element like - for example - cables, board 
20 traces, connectors, packages, sockets, stubs, resistors, inductors or capacitors. 
Even multiple reflections or frequency dependent parameters like skin effector 
dielectric losses can be characterized according to preferred embodiments. 

In testing of modern high-speed electronic devices, e.g. in the frequency range 
of tenth to ten times GHz, the corresponding high-speed signal path e.g. 
25 between an Automated Test Equipment (ATE) and a Device Under Test (DUT) 
becomes relevant influence on the measured signals. In a preferred 
embodiment the device to be modeled is such high-speed signal path, which 



might be connected between the ATE and the DUT. The device to be modeled 
can also be a line drive output of such ATE, an n-port network, etc. 



According to an embodiment, the device to be modeled can be characterized in 
the time domain by at least one of its pulse or step response. An ideal pulse 
5 allows characterizing the device to be modeled in the time domain. Such ideal 
pulse is infinitely short, infinitely high, has a predetermined energy, and is 
known as Dirac pulse. Such Dirac pulse, however, is often not sufficiently 
feasible in reality. Therefore, in a preferred embodiment a substantially ideal 
pulse or the response to a step input signal is used to build the model out of 
1 0 measurement data. 

It is also possible to use the measured external behavior of the device in the 
frequency domain, e.g. using a vector-network-analyzer. For example, the step 
response can be calculated by an inverse Fourier transformation and 
subsequent integration of the resulting pulse response. 

15 BRIEF DESCRIPTION OF THE DRAWINGS 

Other objects and many of the attendant advantages of the present invention 
will be readily appreciated and become better understood by reference to the 
following detailed description when considering in connection with the 
accompanied drawings. Features that are substantially or functionally equal or 
20 similar will be referred to with the same reference signs. 

Fig. 1 shows stimulating incident waves into a 2-port network, 

Fig. 2 shows the SPICE model of the 2-port network of Fig. 1 , 

Fig. 3 shows the model of a voltage adder, 

Fig. 4 shows an example for an approximated step response, 

25 Fig. 5 shows time and voltage values of the sampling points shown in 
Fig. 4, 
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Fig. 6 shows a SPICE model representing a step response, and 

Fig. 7 shows how the four step responses of a 2-port network can be 
measured using a time domain reflection (TDR) oscilloscope. 

MORE DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
5 ACCORDING TO THE INVENTION 

Without limitation of the general inventive concept, in the following the invention 
is described on the example of a 2-port network with ports x and y. As known 
from systems theory such a 2-port network can be fully described by the 
so-called S-parameters in the frequency domain. In an analog manner, the 
10 network can also be described in the time domain by its step response 
according to the following convolution equations: 

e R i(t) = srn(t) * e n '(t) + sr 12 (t) * e l2 '(t) (1) 

e R2 (t) = sr 21 (t) * e M '(t) + sr 22 (t) * e, 2 '(t) (2) 

The term sr xy (t) stands for the response at port x on a stimulation at port y with 
1 5 an ideal step signal. The term ei y '(t) stands for the first derivative of the incident 
wave into port x. The symbol "*" stands for the convolution. 

To put the above equations (1 ) and (2) into a simple recipe: Take each little 
step the incident signal makes, multiply with the step response and delay the 
result by the time when your input step takes place. Summing everything 
20 together then generates the response to the input signal. The responses to 
each individual input signal can be added together according to the principal of 
superposition. 

Fig. 1 shows stimulating incident waves into a 2-port network. The relationship 
between voltages and currents at t he p orts v1 andv2ofthen etwork a re 
25 defined by their characteristic impedances Z1 and Z2. Incident waves can be 
generated by external voltage sources: 
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ens(t) = Zi/(Zsi+Zi)xv S i(t) 



(3) 



e,2s(t) = Z 2 /(Zs2+Z 2 ) x v S2 (t) (4) 
or by reflections at the ports 

e l1R (t) = (Zsi-ZiVCZsi+ZO x e R1 (t) (5) 

5 e, 2R (t) = (Zs2-Z 2 )/(Z S 2+Z 2 ) x e R2 (t) (6) 
According to the principle of superposition, it results from addition: 

e,i(t) = e, 1s (t) + e,iR(t) (7) 

e, 2 (t) = e, 2S (t) + e, 2R (t) (8) 



Thus, from systems theory it is clear how incident waves are generated by 
1 0 transmission line theory and how reflected waves are generated by convolution 
with step responses. This information can be used to construct a 
comprehensive model to be used in the SPICE system. Usually a 50-Ohm 
environment is given and therefore in the following the model is designed with 
50-Ohm characteristic impedance at its ports. Any deviation from that value is 
15 regarded as internal effect and will show up in the step response part of the 
model. 

Fig. 2 shows the SPICE model of the 2-port network of Fig. 1 and Fig. 3 shows 
the model of a voltage adder. 

Calculating the incident waves from the schematic result into: 



20 e M (t) = (Z, x v S i(t) + (Z S i - ZO x e R1 (t))/(Z S i + Zi) (9) 

e, 2 (t) = (Z 2 x v S2 (t) + (Z S2 - Z 2 ) x e R2 (t))/(Z S2 + Z 2 ) (10) 

e R i(t) = srn(t) * en'(t) + sr 12 (t) * e l2 '(t) (1 1 ) 

e R2 (t) = sr 21 (t) * en'(t) + sr 22 (t) * e, 2 '(t) (12) 
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Using simple mathematics, it can be shown that these equations match with the 
calculated equations mentioned above. 

Any arbitrary step response can be approximated with arbitrary accuracy by a 
section-wise linear curve like in the example shown in the top graph of Fig. 4. 
5 The corresponding sampling points for the approximation in this example have 
the time and voltage values listed in Fig. 5. In order to construct the step 
response, the incoming step s(t) is stripped down to the pulses Pi(t) having their 
transitions at the sampling points used for approximation. Those pulses are 
used to control current sources with a trans-conductance representing the 
10 slope in the corresponding section of the step response like depicted in Fig. 6. 

The delays TDj of the transmission lines T s (i = 1 ¥ 2 8) represent the lengths 

of the construction pulses pj(t): 

TDjMj-tM (13) 

The differential input signals on the voltage controlled current sources Gj look 
1 5 like the construction pulses pi in Fig.4. The trans-conductances GTj have to be: 

GTi = (Vi-V M )xC1/(ti-ti) (14) 

The resistors R1 and R2 define the DC path for the output voltage. They can 
be calculated from the end value of the step response: 

R2/(R1 + R2) = sr(oo)/s(oo) = sr(oo)/1V (15) 

20 The time constant of the RC element 

tr C = C1 x(R1 II R2) (16) 

has to be large enough so that the exponential loading curve does not interfere 
with the shape of the step response. A good value would be if the time constant 
were 1 000 times larger that the time range between the first and the last 
25 sampling point. This limits the resulting voltage error of the simulated step 
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response to below 0.1%. 

With the SPICE model shown in Fig. 6 a circuit is created that has the step 
response needed for the device to be modeled. As the model is linear and time 
invariant, its response to any arbitrary input signal is equal to the response of 
5 the signal path to be characterized in the described embodiment. 

Fig. 7 shows how the four step responses of a 2-port network can be measured 
using a time domain reflection (TDR) oscilloscope. The internal TDR source is 
dedicated to the first channel CH1 . It generates a steep edge with a rise time in 
the order of 40ps. This can be taken as an approximate ideal step scaled with a 
10 factor. The delay time relative to the internally generated oscilloscope trigger 
can be observed as the point in time on the scope graph of channel 1 where 
the edge starts. 

Setup A places a short circuit at the end of cable 1 generating a reflection 
coefficient of -1. Channel 1 of the oscilloscope sees a pulse with a length that 
1 5 equals the round trip delay of cable 1 . The delay of single trip of cable 1 can be 
easily determined by placing markers at the 50% points of the pulse edges and 
taking half of the time distance between the two markers. 

In Setup B the two channels capture the same signal but on channel 2 (CH2) 
the combined delay of cable 1 and cable 2 takes effect. The delay of cable 2 
20 can be determined by placing markers at the 50% points of the two channels, 
reading out the time distance and subtracting the delay time of cable 1 . 

The Setups A and B can be regarded as a calibration procedure. 

In Setup C and Setup D the reflections at the stimulated port can be observed 
as superposition to the stimulating time domain reflection pulse on channel 1 of 
25 the oscilloscope with an additional roundtrip cable 1 delay. The transmitted 
signals can be observed on channel 2 with an additional delay of cable 2. 

The four step responses can now be derived very easily from the measurement 
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and from using equations representing the Setup A to Setup D described 
above. 

Although the edge of a TDR edge is usually very fast, it certainly is not a 
perfect step. In order to eliminate the remaining distortions due to finite slew 
5 rate, possibly over- and undershoot and settling phase to final value, the ideal 
step response can be calculated from the measured stimulus and response by 
the method described below. 

At first, the captured waveforms together with a shape of an ideal step starting 
at the same time as the stimulus have to be converted into a periodic shape by 
10 extending the measured date set with the same edge that has been mirrored at 
the average value of the first and the last sampling point. 

The resulting signals can then be converted into the frequency domain via Fast 
Fourier Transformation (FFT). The transfer function results from the division 
from the output signal by the input signal. Multiplication of transfer function and 
15 the ideal step in the frequency domain gives the spectrum of the ideal step 
response that can be converted back into time domain via inverse FFT. 

As the output voltage of the model shown in Fig. 2 is generated by integrating 
loading currents, the voltage error in the transient analysis strongly depends on 
the integration method and the time steps of the SPICE simulator. As a rule of 

20 thumb, the time step should be in the same order of magnitude as the delay 
time of the transmission lines in the step response model. Depending of which 
simulator is used, there may be a lot more options controlling simulation 
accuracy. The best way to test those parameters is to stimulate the model with 
a single edge and verify if the final value is close enough to the theoretical end 

25 voltage. 

Differing from the comprehensive model shown in Fig. 2 there are applications 
where not all step responses are of interest and for such applications the model 
can be drastically simplified. In many cases only the transmission of a signal 
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from port 1 to port 2 is necessary. This is the case for example for an 
environment with 50-Ohm terminations at both ends or for signal paths with 
only small reflections. In such cases, the responses sr1 1 , sr21 , and sr22 can 
be neglected and deleted from the model. The measurement effort and the 
5 model complexity are strongly reduced. 

Another example for the device to be modeled is a line driver controlled by a 
digital input signal. In this case, there is no response at the control input to 
anything happening at the output. Therefore s12 and s1 1 are not necessary. 

Similar to S-parameters in the frequency domain, step responses can be used 
1 0 to characterize n-port networks with n > 2. For example, coupled transmission 
lines for differential signaling would be an application for it. 

In summary, the present invention introduces a new method of modeling linear 
and time invariant devices like signal paths or n-port networks. The principle 
has been derived from basic transmission line and systems theory and is based 
1 5 on step response measurement with a TDR oscilloscope. It has been validated 
by the means of a SPICE simulation showing the correlation between an 
arbitrarily chosen signal path model and the corresponding step response 
model. 
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